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ABSTRACT 


A previous statistical survey of data from the HPCE experiment on the AMPTE/CCE 
satellite established probability distributions for trapped ions and electrons. An extension 
of this survey for ions at 240 and 442 eV and for electrons at 340 and 770 eV confirmed 
these distributions. A further detailed analysis of the electron data firom 13 individual 


data collection days also showed the trapped electron distributions to be concentrated in 
the dawn to noon region, centered at L = 7. These trapped electron distributions can be 
described as a bi-Maxwellian distribution function and be characterized reaeuuably by khe 
criteria that the flux has to exceed 5x10* (cm* s sr)'\ the distribution has to be within 10* 
of the magnetic equator, the ratio of the perpendicular temperature to the parallel 
temperature is greater than 3 and that the anisotropy is greater than 2.0 for 150 eV 
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I. ZMTRODDCTION 


Equatorially trapped plasmas are ion and electron 
distributions trapped within a few degrees of the Earth's 
magnetic equator. Equatorially trapped plasmas can be 
described by a bi-Maxwellian distribution function. These 
trapped plasma distributions were defined by the initial 
observations by Olsen (1981). Ion and electron distributions 
with highly anisotropic pitch angle distributions, peaked at 
90 degrees pitch angle, were observed at energies from a few 
eV to hundreds of eV, near geosynchronous orbit. These trapped 
distributions are of interest as indications of basic wave 
particle interactions, and as an intermediate process in 
plasmasphere filling. 

The energy - pitch distribution indicates the wave 
particle interaction aspect of these plasma distributions. 
They are indicative of perpendicular acceleration (Tp„p > Tp,r) , 
and quasi-linear diffusion (flat diffusion at low energy). 
Though not yet proven, there are indications of a 
correspondence between equatorially trapped plasmas and 
Bernstein mode waves (equatorial noise) and electron cyclotron 
harmonics (Gurnett, 1976 and Kurth et al., 1979). 

The plasmasphere filling role is indicated by the 
correspondence between the plasmapause region and the location 
of equatorially trapped ion distributions (Horwitz et aJ., 
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1981). The variations in pitch angle structure with latitude 
also suggests the role (Olsen et al., 1987). 

There have been previous surveys of equatorially trapped 
plasmas. Olsen et al., (1987) surveyed DE 1/RIMS (ion) data 
for 0 - 100 eV. Sagawa et al. (1987) surveyed the DE 1/EICS 
(ion) data for 0-1 keV. Both surveys of equatorially trapped 
plasmas were limited in altitude by the DE 1 orbit, which had 
apogee at L = 4.7. Both of these surveys also lacked 
complementary electron data. 

Braccio (1991) surveyed AMPTE/CCE data for both ions (30 - 
150 eV) and electrons (150 eV) out to 8.8 Earth radii. 

The first purpose of this thesis is to extend the survey done 
by Braccio by conducting a statistical survey of the next two 
higher energy channels for both the ions and electrons (i.e. 
240 and 442 eV for the ions and 340 and 770 eV for the 
electrons.) The survey will be used to extend the findings in 
his paper and to attempt to answer questions he raised, such 
as whether the higher energy ions inhabit the region of 
decreased probability from 1200 to 1400 local time. The 
electron survey will ensure that the trapped distribution is 
not being under-evaluated and that the location and shape of 
the trapped distribution found by Braccio is accurate. 

The second purpose of this thesis is to conduct a detailed 
analysis of 13 individual days to examine the nature of 
equatorially trapped distributions. The trapped distributions 
will be described as bi-Maxwellian distributions, in order to 
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quantify a trapped distribution. In addition, this thesis will 
be used to better define the electron distributions, and the 
location of the trapped electron distributions. 
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ZI. BACXOSODND 


A. THE PLASKASPHEKE 

A magnetosphere is the region around a magnetized 
planetary body in which that body’s magnetic field plays the 
dominant role in defining the behavior of charged particles. 
It's outer boundary, the magnetopause, occurs where the solar 
wind, and the magnetic field in the solar wind, becomes 
dominant. This boundary occurs on the sunward side of the 
Earth at approximately 10 Earth radii (roughly 63,750 km). The 
location of this boundary is determined by a balance between 
the pressure exerted by the solar wind and the obstacle formed 
by the Earth's magnetic field. During active times, the 
magnetopause has been observed as close as 5 geocentric Earth 
radii (Re). The inner boundary of the magnetosphere occurs at 
the top of the ionosphere. This boundary can be taken as 
occurring at an altitude of 1000 km or 1.16 R^. (Parks, p.7) 

As can be seen in Figure 1, the Earth's magnetosphere is 
highly asymmetric. While the sunward boundary is located at 
approximately 10 Re, the Earth's magnetic tail has been seen 
to extend beyond 200 Re on the nightside. The length and shape 
of the magnetic tail again depends on the interaction between 
the geomagnetic field and the solar wind. (Parks, pp. 7-8) 
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For our purposes, the major components of the 
magnetosphere are the plasmasphere, the plasmasheet, and the 
plasmapause. The plasmasphere is the region of the 
magnetosphere that is closest to the Earth. It begins just 
above the ionosphere, at low to mid-latitudes. The 
plasmasphere extends in altitude to between 3 and 5 in the 
equatorial plane, and between +- 60“ magnetic latitude just 
above the ionosphere. The plasmasphere corotates with the 
Earth and particles in this region are affected by the Earth's 
corotational Electric field. (Parks, pp. 11 and 73) 

This region contains plasma, ionized atoms and electrons, 
with densities of 10^ - 10^ cm‘^. Characteristic ion and 
electron energies are on the order of l eV at 4.5 Rg. The 
density of the plasmasphere decreases with altitude. In 
general, the density in this region experiences a gradual drop 
proportional to the fourth power of the Mcllwain L parameter 
(a measure of altitude based on the magnetic field lines that 
will be discussed later). This is illustrated in Figure 2. 
(Chappell et aJ., 1970) 

At approximately 3 to 5 Earth radii, again depending on 
the magnetic activity history, the plasmapause is encountered. 
This is a transition region for the plasma in which plasma 
energies sharply increase (Parks, pp. 231 and 502) and 
densities drop, generally very sharply, that is used to define 
the inner boundary of the plasmapause. Figure 3.(Harris et 
al., 1970). 
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Composite of the first four H* ion concentration profiles for the inbound and out¬ 
bound passes of OGO 5. 


Figure 3. Plasmapause Magnetic Activity Dependence 
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These aspects of the plasmasphere can be further 
illu'=trated using (relatively) more modern data from ISEE 1 
total electron density measurements obtained from observations 
of plasma waves. Figure 4a shows density versus L, with a 
solid line at 100 x (L/4.5)'^ superimposed. The plasmapause is 
at L = 4.8, 1700 local time. The plasma density outside the 
plasmapause continues to drop as L~*. This characteristic of 
the plasma density profile can more easily be seen if the data 
are normalized by V, as in Figure 4b. (Olsen, 1992) 

The plasmasphere density is dependent on the magnetic 
activity. A large magnetic storm can effectively push the 
plasmapause in to less than 3 R^. Figure 3 shows the effects 
of magnetic activity on the density and location of the 
plasmapause. Magnetic intensity increases from a low in the 
upper left panel in the figure to a maximum in the lower right 
hand panel. (Harris et ai., 1970) 

The storm-time electric field strips away the plasma at 
higher altitudes, as the plasma is convected to the 
magnetopause. This region is then refilled from the ionosphere 
after the storm-time field relaxes. The process, termed polar 
wind, is driven by ambipolar diffusion after the electric 
field relaxes back to a steady state value of approximately 1 
mV/m. This diffusion process calls for electrons to leave the 
upper ionosphere, probably driven by photoemission, and move 
along the geomagnetic field lines. The resulting ambipolar 
electric field, caused by the displacement between the 
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electrons and O* in the upper atmosphere, causes lighter ions, 
such as and He% to be dragged up the field lines after the 
electrons. This 'polar wind* results in a refilling rate of 1 
to 10 ions/cm^ per da>. (Horvitz, 1983) 

On the nightside the plasmasphere is bounded by the 
plasmasheet, a region of low density, hot plasma (with 
densities on the order of 1 cm'^ and characteristic energies 
of 1-10 keV) . The corresponding plasmapause for this region is 
very distinct and the transition from plasmapause to 
plasmasheet takes place rather rapidly. This is not the case 
for the region that extends from just before dawn until just 
after dusk local time, on the dayside. (Parks, pp. 231 and 
502) 

On the sunward side of the Earth, the plasmapause is a 
region that can be as much as 1 Re in width. Additionally, 
there is usually no sharp distinction corresponding to its 
inner and outer boundaries during this local time period 
(Parks, p.231). Therefore, it is usually a matter of judgement 
as to which region you are studying. 

The region between the dayside plasmapause and the 
magnetopause is also ill defined. It is not clear whether the 
plasmasheet encircles the Earth and occupies this region. 
While there is no known reason why this should not be the 
case, the plasma observed in this region does not display the 
characteristics of that which is found in the nightside 
plasmasheet. This has led to questions concerning the plasma 
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filling mechanism for this region as well as to questions of 
where the plasma in this region comes from. 

In the dusk region there is an additional asymmetry as 
seen in Figure 5. This dusk bulge is the result of interaction 
between the corotational electric field of the plasmasphere 
and the cross tail electric field induced by the solar wind. 
The corotational field is the result of the charged particles 
rotating with the Earth while trapped in its geomagnetic field 
and is directed radially inward toward the Earth. The cross 
tail electric field is induced by the solar wind's interaction 
with the Earth's magnetic field. This cross tail field is in 
the dawn-dusk direction in the equatorial plane of the 
magnetotail. The sum of these two electric fields results in 
a series of equipotential contours which mirror the dusk 
bulge, Figure 6. (Parks, pp. 231-236) 

B. THEORY 

1. Plasma Definition 

A plasma is a collection of discrete ionized and 
neutral particles, which has overall electrical neutrality. 
The physical dimensions of the plasma must be large in 
comparison with a characteristic length lo called the Debye 
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Equipotentlal contours for the magnetospheric electric field in the 
equatorial plane. Upper left: first-order approximation lor tfie convec t ion electric 
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Figure 6. Magnetosphere’s Electric and Magnetic Fields 
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length. The total number of charges in a sphere with a radius 
Xd must be much greater than 1. 

2 . Debye Shielding 

The electrostatic potential of a point charge g in a 
vacuum is given by: 

(1) —2— Volts 

where f is the electrostatic potential, e,, is the permittivity 
constant and r is the distance from the point charge (Halliday 
and Resnick, 1988) . If the charge is immersed in a plasma, a 
positive charge will attract electrons while repelling ions, 
and similarly, a negative charge will attract ions and repel 
electrons. The potential then becomes: 

(2) ,|r=._S_exp"^ 

where k is the Boltzman constant, T, is the electron 
temperature, which is a measure of the average kinetic energy, 
n is the equilibrium density of the plasma, and e is the 
charge of the electrons. This has the effect of screening out 
electric potentials in a plasma. The electron temperature is 
used in the definition of because the electrons are more 
mobile than the ions and do most of the shielding by creating 
a surplus or deficit of negative charge. (Parks,1991) 
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3. Plasma Parameter 


In order for a collection of ionized particles to be 
considered a collisionless plasma, three conditions must be 
satisfied. The Debye length must be much less than the 
dimensions of the plasma. The number of particles in a Debye 
sphere defined as: 


(3) 






must be much greater than 1, that is, there must be a large 
enough number of particles for Debye shielding to be 
statistically valid. Finally the frequency of collisions of 
particles must be low. The plasma parameter g is defined as: 


(4) 

and for a collisionless plasma g « 1. (Parks, 1991 and Chen, 
1983) 

4. Motion in a Uniform Magnetic Field 

A collisionless plasma will behave as a collection of 
individual particles. The individual charged particles will 
move in trajectories determined by the applied electric and 
magnetic fields. For space applications the fields resulting 
from the particle motion are often small and may be neglected. 
This is particularly true for the magnetic field, less so for 
the electric field. The force acting on a charged particle 
moving in a combined E and B field is given by the Lorentz equation: 
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(5) 




where £ is the Lorentz force, and v is the particle velocity. 

For the case where £ - 0, and the magnetic field is 
uniform, a charged particle will execute simple cyclotron 
gyration with a frequency: 


( 6 ) 


“c= 


igi-B 

m 


and radius: 


(7) 






_ ^perp _ invsintt 
\glB ~ \g\B 


where m and v are the mass and velocity of the charged 
particle and the pitch angle a is the angle between the 
velocity vector of the charged particle and the magnetic 
field. The velocity parallel to the magnetic field line, Vp„ 
= vcoso, is not affected by the magnetic field. This motion 
describes a circular orbit about a guiding center which is 
travelling along the magnetic field line with velocity Vp„. 
The trajectory of this particle is a helix with it's axis 
parallel to the field line. (Chen, 1983) 

5. Magnetic Mirror 

The Earth's magnetic field lines converge at both the 
north and south magnetic poles, and the field strength 
increases with altitude. Because of this, there is an 
additional force that acts upon the charged particle. This 
force can be expressed as: 
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( 8 ) 


£=-|iVB 



where m is the magnetic moment. Because the gradient in the 
magnetic field is parallel to the direction of the field line 
it can be seen that this force is directed along the field 
line. Since the force is directed parallel to the particles 
parallel velocity component, it will obviously affect the 
particle's velocity. From Lentz's law, it can be shown that n 
is invariant. (Parks, pp. 89-90) 

Since n is invariant, Vperp must increase as B is 
increased. For this to happen Vp^j. must decrease since v^perp = 
v^ - v^par (from conservation of energy) . Therefore, given a 
large enough B, there will come a point where Vpar == 0. At this 
point Vpajp will equal v, and the particle will mirror back 
along the field line. Figure 7. (Gladstone, 1967). (Notice 
that the gyroradius also gets smaller as the particle 
approaches the mirror point as a result of its Vp^^p 
dependence.) 

For a particle that mirrors, equation 8 then leads to: 

2 2 

jgj ^^perpo _ ^^perpm 

QBo gB„ 

where subscript o refers to values at the equator and m to 
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those at the mirror point. Rearranging equation 9 gives: 


( 10 ) 


^o_ ^p9xpo ^ ^pezpo 

v| 


defining the pitch angle of a particle, a, to be the angle 
between velocity vector of the partide and the magnetic field 
line gives Vp^^p = v sin a. Plugging this into equation 10 
gives: 

(11) :^=sina| 

which states that all paz^icles with a pitch angle Sq will 
mirror at the location defined by B = Bi,. Particles with a > 
tto will mirror at lower altitudes. (Parks, pp. 111-112) 

The magnitude of the Earth's magnetic field is given 
according to: 


( 12 ) 


S=S 


OS 


V4-3cos^X 

L^cos®A 


where B^, is the magnitude of the Earth's magnetic field on the 
Earth's surface at the magnetic equator, X is the magnetic 
latitude, and L is the Mcllwain L parameter (Parks, p. 54). 
The Mcllwain L parameter is variable, given in units of Earth 
radii, used to label magnetic field lines with relation to 
where they cross the plane of the magnetic equator. Its value 
is given by: 
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(13) 


cos^A 

where r is the distance from the Earth center, in Re, to the 
field line at the magnetic equator. (Par)cs, p. 115) 

Sxibstituting equation 12 into equation 11 gives: 


(14) 


sin^aa= 


cos*Aa 

^-Scos^Aj, 


Therefore, by defining an equatorially trapped plasma to 
mirror at a magnetic latitude of ± lOo or less, this requires 
that a charged particle have a pitch angle greater than 69°, 
at the equator, in order to be equatorially trapped. Figure 8. 
It is these trapped particles that will be investigated in 
this paper. 

6. Statistical Distribution 

When dealing with a system which is composed of a very 
large number of individual particles it becomes impractical to 
solve the equations of motion for the system. Instead the 
particles may be treated statistically through the use of the 
distribution function. Classically the distribution function 
gives the probability distribution of the values of the 
coordinates and momenta of the particles. The density of 
particles in coordinate space is obtained by integrating the 
distribution function over the momentum or velocity space. 
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Pitch Angle vs Magnetic Latitude 



Figure 8. Pitch Angle verses Mirror Latitude 
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That is: 


(15) ■n=j*.f (X/if) 

where n is the density, and f(£,y) is the distribution as a 
function of position and velocity integrated over three 
dimensional velocity space. 

The two distribution functions of interest in this 
thesis are the bi-Maxwellian, and the Maxwellian or isotropic 
distribution function, given by: 



where T p,jp and Tpa^. are the characteristic temperatures in the 
perpendicular and parallel directions with respect to the 
magnetic field line. 

C. PREVIOUS OBSERVATIONS 

Thermal plasma pitch angle distributions seem to have been 
first studied by Horwitz and Chappell (1979) and Comfort and 
Horwitz (1980) . These authors used electrostatic analyzer data 
to study ion pitch angle distributions at geosynchronous 
orbit, using ATS 6 data taken in 1974. The surveys dealt with 
data taken at 10.50 off the magnetic equator. 
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Comfort and Horwitz (1980) observed two important aspects 
of ion pancake distributions (peak flux near 90o pitch angle) . 
The first was the occurrence probability for the pancake 
component of the ion distribution was local time and energy 
dependent. The highest probability of occurrence occurred in 
the lowest energy channel (20 - 40 eV) studied and for local 
times between 1400 and 1800. Pancake distributions wez«3 seen 
42% of the time in this sector for ions of that energy. 

The second was that Comfort and Horwitz observed that 
field aligned ions and ions with 90 o pitch angle seem to be 
anti-correlated at 1700. Figure 9 shows that there is a 
decrease in the occurrence probability of field aligned ions 
when there is a peak in the pancake occurrence probability. 

Horwitz et aJ. (1981) studied pancake distributions in low 
energy (^ 100 eV) ion data obtained from the ISEE 1 mass 
spectrometer. The H^ distributions were often found in the 
vicinity of the plasmapause,Figure 10, and usually just inside 
the plasmapause. Horwitz et al. also observed that the pancake 
distribution was often seen in the presence of colder, 
isotropic plasma. 

Olsen (1981) observed a thermal plasma population, trapped 
within a few degrees of the magnetic equator, using 
electrostatic analyzer data from the SCATHA satellite. Figure 
11 shows the ion count rate, for various ion energies, as a 
function of pitch angle. The data for this plot was taken at 
the equator at approximately 1000 local time and 5.5 Rg. This 
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Tottil percent occurrence frequencies 
for ion pitch angle distributions having 
deilgnated components, as functions of local 
time* 


Figxire 9. Field Aligned and Pancake Trapped Ion Distributions 













MTCH A .OlC (®» 

Ion pitch angle distributions st the e«|ustor for dsf 179 of 1979. Data are 
plotted froM the FIX ahd HI detectors at 11, dt, 193, 523 , and 900 eV. Fluxes are 
sealed by increasing factors to keep the* frote overlapping. Maxlaua values,^with^ 
Increasing energy are 1150 c/s, 7300 c/s, 9200 c/s, 550 c/s, and 590 e/s. The 0*-180* 
range corresponds to looking sunward, with 0* corresponding to looking south. 


Figure 11. 


ion Pitch Angle Distribution 
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figure clearly shows a trapped distribution, centered at 90© 
and 2700 pitch angle, for ions of energies 11 to 103 eV and, 
to a lesser extent, for those at 523 eV. The 900 eV ions do 
not show evidence of a trapped distribution. This figure also 
shows a well defined loss cone for the three highest energies. 

Olsen observed a like distribution in the electron data. 
Figure 12. A source cone, centered at Qo and ISO® pitch 
angles, was seen in the 41 eV electron flux concurrent with 
the trapped distribution at higher energies. This led to 
speculation that the field aligned particles were the 
(ionospheric) plasma source, and these particles were 
subsequently heated in the traverse direction. Note that the 
count rates in the last two figures are scaled differently for 
different energy levels in order to facilitate presentation of 
the data. Figure 13, from Scott (1991), shows a plot of count 
rate versus energy (in eV). The trapped electron distribution 
is seen to exist in the 50 to 1000 eV range, corresponding to 
temperatures of 100 - 200 eV and densities of 1 - 10 cm‘^. The 
trapped ions show a peak in the 20 to 200 eV range. This 
corresponds to temperatures of 20 to 50 eV and densities of 1 

- 10 cm'\ Sagawa et al. (1987) observed a local time 
dependence in the location of the trapped ions in data from 
the Dynamics Explorer (DE) 1 satellite. Sagawa et al. 
additionally saw that the trapped ions were composed primarily 
of H* ions and that these were in the lowest energy bin (0.01 

- 1 keV) of the DE 1/EICS summary plots. They reported that 
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Electron pitch angle distributions at 
the equator for day 179 of 1979. Data are plot¬ 
ted from the HI detector at 41, 523, and 4730 eV. 
Pitch angle conventions are as in Figure 7. 


Figure 12. Electron Pitch Angle Distribution 
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Ion and electron count rates as a func¬ 
tion of energy from the LO detector near 90* 
pitch angle for day 179 of 1979. The electron 
count rate has been scaled by a factor of 2. The 
count rates from the FIX detector dwells were se¬ 
lected at their maxima (90* pitch angle). The 
difference between the LO and FIX ion data re¬ 
flects degradation of the spiraltrons for the LO 
ion detector. The LO count rate curve has been 
traced and moved up to overlap the FIX detector 
data (about a factor of 2). The peak in ion 
count rate at 700 eV is a local maximum in the 
distribution function as well (see Figure 10). 


Figure 13. Energy Relations of Trapped Plasmas 
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the Hcllwain L value was higher, for the peak ion occurrence 
probeibility, in the local noon and dusk sectors than it was 
near local midnight, Figure 14. Olsen et al. (1987) also saw 
this in their statistical survey. Olsen et al. noted that the 
latitudinal extent of the high probability region is local 
time dependant, ranging from ± 30 o in the early afternoon 
region to ± 10° in the early dawn region. 

Olsen et al. (1987) observed, from data collected by DE 1, 
that the trapped ion distribution was composed primarily of 
H*, but that He* was seen to have a trapped component, having 
10% the density of the trapped H*, approximately 40% of the 
time. In one case, trapped O* was seen with a relative density 
of 0.1% that of H*. Additionally, the trapped distribution was 
observed to be very localized about the equator. This is seen 
in the fact that the ions change from a field aligned 
distribution to a trapped distribution and then back very 
quickly as the satellite traverses the equatorial region. 
Figure 15 illustrates this aspect of the evolution in the 
pitch angle distributions. 

Figures 15a, 15b, and 15c show plots of flux verses pitch 
angle for the magnetic latitudes of -7.9°, - 1 . 9 o 
(approximately), and 3.6° respectively. In this case, the He* 
ions mirror the H* ions, although at about 3.5% of its flux. 
Figures 15d, 15e, and 15f show the distribution functions for 
H* in these time periods. Notice the drop in density and the 
increase in temperature as the satellite enters the equatorial 
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Occurrence probability of low>enerfy (0.01-1 keV) H'*' 
pancake distribution with peak ion Bux above 10^(cm^t sr)~^, 
within 5^' of the magnetic equator for active times (A'p^3-) as a 
function of MLT and L shell. L shell bin size AL=0.5. 


Figure 14. Trapped Ion L verses Local Time Dependence 
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region. (Olsen et al., 1987) 

Klumpar et al. (1987) found examples of eguatorially 
trapped plasma in the data from the ANPTE/CCE satellite. The 
trapped ion distribution was found near the plasmapause 
interface. The temperature of these ions were found to be on 
the order of 30 - 50 eV. Like Olsen (1981), they also observed 
that the angular distributions of the trapped ion distribution 
beccune narrower for increasing ion energies. 

Braccio (1991) extended the survey of the AMPTE/CCE data 
by surveying ions in the 30 - 150 eV energy channel and 
electrons in the 150 eV energy channel out to 8.8 Earth radii. 
In this survey Braccio concluded that the location of 
eguatorially trapped plasmas is species dependant. The ions 
and electrons show a different local time dependence in the 
location of their occurrence probability peak. The following 
summarizes his conclusions: 

Electrons 

• display a uniform high probability distribution centered 
at 0900 local time and a L value of 6 

• display a weak, if any, L dependence on local time 

• trapped electrons begin to be seen at dawn leading to 
speculation that their existence is dependent on 
photoelectron emission from the Earth's ionosphere 

• the shape of the distribution is basically conical, 
however, it drops off more rapidly than it decreases, with 
respect to local time 
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Ions 


* distribution shows a strong L dependence on local time 

• high trapped ion probability region begins at local dawn, 
for L approximately 4, and rises to a maximum at between 
1400 and 1500 local time with an L value of 8 

• the distribution drops off quickly in altitude after 
peaking as local time increases there is a region of 
decreased probability in the afternoon sector for his data 
that he suggested may be inhabited by higher energy ions. 
Part of this thesis will investigate this. 

* the overall shape of the high trapped ion prob 2 d}ility 
region mirrors that of the location of the plasmapause 

' trapped electrons seem to be excluded from regions of high 
trapped ion probabilities and vice versa. 

This thesis will further extend the look at the AMPTE data 
by conducting a statistical survey of ions in the 240 eV and 
442 eV energy channels and electrons in the 340 eV and 770 eV 
energy channels. 

There is not always a clear criteria used to define an 
equatorially trapped plasma. Intuitively, the definition seems 
to be that the bulk of the plasma is confined to (mirror 
within) 10° of the magnetic equator. One aspect of the 
development which follows is an attempt to quantify this 
concept, and relate the various definitions, tests, and survey 
techniques. Braccio (1991) defined, for the purpose of his 
study, an equatorially trapped plasma as having a minimum 
(threshold) flux in the 80® to 90° pitch angle bin of 10® (cm^ 
s sr)‘^ for ions and 5 x 10® (cm^ s sr)‘^ for electrons with an 
anisotropy greater than 1.5 (The anisotropy is defined as the 
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ratio of the fluxes in the 80o to 90o pitch angle bin to those 
in the 60o to 70o pitch angle bin.) in the same period. His 
survey generally was restricted to data taken within 10° of 
the magnetic equator. For the purpose of the statistical 
survey in this paper, the anisotropy criterion of 1.5 remains 
the scune for both ions and electrons and the minimiim flux 
criterion for a equatorially trapped plasma will be as 
follows: 

Ions 

• 240 eV: 10‘ (cm^ s sr)‘^ 

• 442 eV: 10^ (cm^ s sr)"^ 

Electrons 

• 340 eV and 770: 5x10* (cm^ s sr)’^ 

The lower flux values for the ions partly accoxuits for the 
lower flux of ions at the higher energy channels. 

Another criterion which could be used to characterize a 
trapped distribution is the ratio of the plasma's 
characteristic temperatures in the perpendicular and parallel 
directions with respect to the magnetic field line; Tp^^p and 
Tpari. An equatorially trapped distribution may be described by 
a bi-Haxwellian distribution function with these 
characteristic temperatures (Scott 1991). For the purpose of 
the case studies of this thesis a ratio of Tp^^p to Tpa^i equal 
to three will be used as a lower threshold required for an 
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eguatorially trapped distribution to occur. The relationship 
berveen the above two definitions; flux ratios and temperature 
ratios, is considered next. 

Applying the threshold of a temperature ratio of three 
determines the anisotropy in flux that must be met for an 
electron distribution to be considered trapped. Figure 16 
shows plots of the log of the flvix versus pitch angle for 150 
eV and 340 eV. In each the perpendicular temperature was set 
to 75 eV and the parallel to 25 eV; a temperature ratio of 
three. These temperatures were selected as representative of 
values seen in the data analysed. Scott (1991) got slightly 
larger ratios. The density was set arbitrarily to 1.2 cm'* as 
a reasonable value. The density does not affect the 
illustration. 

In both Figure 16a and 16b there are shown two curves. The 
uppermost is the plot of the bi-Maxwellian distribution 
function with the above Tp,tp and Tp.^i values at the magnetic 
equator. The lower of the two curves is a plot of the same 
function with the same temperatures, but mapped to a magnetic 
latitude of 10'’. The peak flux of this curve is equal to the 
magnitude of the flux at the magnetic equator at a pitch angle 
of 69°. This is because a pitch angle of 69° at the equator 
maps to a pitch angle of 90° at a magnetic latitude of 10°. 
With our previous criterion for a trapped distribution being 
one that is limited to within 10° of the magnetic equator this 
establishes a lower limit on the pitch angles analyzed at the 
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equator at 69° and a relevant pitch angle range of 69° to 90° 
for the analysis of a bi-Maxwellian distribution. 

Taking the ratio of the top curve to the bottom curve in 
Figures 16a and 16b a differential anisotropy was derived. The 
results of this calculation was 1.67 for 150 eV and 3.2 for 
340 eV. The survey work done here, and those by Braccio (1991) 
use an integral anisotropy calculated by taking the ratios of 
the flux in the 80° to 90° pitch angle bin to the flux in the 
60° to 70° pitch angle bin. Note that the survey files were 
divided into pitch angle bins of 10° in size. Using this 
definition for the anisotropy we take the ratio of the 
integral of the distribution function over pitch angles 80° to 
90° to the integral of the distribution function over pitch 
angles 60° to 70°. The results of this is an anisotropy of 
1.953 for 150 eV and 4.421 for 340 eV. For the purposes of the 
detailed analysis case studies in this thesis the rounded 
values of 2.0 and 4.4 will be used for 150 e'^ and 340 eV , 
respectively, as applied to a definition of being bi- 
Maxwellian and hence an equatorially trapped electron 
distribution. This will allow for comparison between this work 
and previous survey techniques. Note that the integral ratio 
calculated here is higher than that used in the initial 
survey. 
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Figure 16. Log of Flux verses Pitch Angle - 150 eV and 340 eV 
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D. THE EMPTE/CCE SATELLITE 

The Active Magnetospheric Particle Tracer Explorers 
(AHPTE) mission consists of three satellites that were 
launched on August 16, 1984. The purpose of this mission 

(Acuna et al., 1985) was to: 

1. investigate the transfer of mass and energy from the 
solar wind to the magnetosphere and to study its further 
transport and energization within the magnetosphere. 

2. to study the interaction between artificially injected 
and natural space plasma. 

3. to establish the elemental and charge composition and 
dyneunics of the charged population in the magnetosphere 
over a broad energy range. 

Two of the satellites, the Ion Release Module (IRM) and 
the United Kingdom Sub-satellite (UKS), were concerned 
primarily with the introduction of artificially injected ions 
into the magnetosphere and will not be discussed further. The 
third satellite was the Charged Composition Explorer (CCE), 
Figure 17. The purpose of this satellite was to measure the 
particle distribution of the naturally occurring plasma, with 
respect to species, energy, and pitch angle, as well as to 
measure the artificially released ions from IRM. (Dassoulas et 
al., 1985) The CCE was placed in an elliptical orbit around 
the Earth with a period of 15.66 hours and an inclination to 
the Earth's equatorial plane of 4 . 8 o. it had an altitude at 
perigee of 1108 km and at apogee of 49,684 km (roughly 1.2 and 
8.8 R£ respectively) . It was spin stabilized with a spin rate 
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of 10.25 rev/min (Dassoulas et al., 1985). Data began to be 
collected by this satellite on August 26, 1984. 

The payload of the CCE, Figure 18, consisted of five 
experiments; 1} the Hot Plasma Composition Experiment (HPCE), 
2) the Charge Energy Mass Spectrometer (CHEM) , 3) the Medium 
Energy Particle Analyzer, 4) the Magnetometer, and 5) the 
Plasma Wave Experiment (Dassoulas et al., 1985). This paper 
concerns itself with data from the HPCE (and, indirectly, the 
Magnetometer). 

E. THE HOT PLASMA COMPOSITZOH EXPERIMENT (HPCE) 

The HPCE consists of the Ion-Mass Spectrometer and the 
Electron Background Environment Monitor (EBEM). The ion-mass 
spectrometer. Figure 19, provides mass per charge ion- 
composition measurements from very low energies (corresponding 
to the spacecraft potential) to approximately 17 keV. The ions 
enter the detector through a collimator which limits both 
azimuthal and elevation angles of acceptance. The azimuthal 
limits are at ±5.5® while the elevation acceptance angle 
ranges from approximately ±25© for ions at the spacecraft 
potential to ±7.5© for those at 17 keV. The ions are sent 
through a retarding potential analyzer (RPA) and then 
accelerated through a -2960 V potential. (Shelley et al., 
1985) 
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Figure 19. The HPCE Ion-Mass Spectrometer 











The ions then pass through the object slit and into the 
cylindrical electrostatic energy analyzer. The electrostatic 
energy analyzer is programmadsle in 32 energy per charge steps 
from 3 to 20 keV/e. The central portion of the ion flux then 
enters the mass analyzer through a second slit, with a portion 
of the spectrum measured by "energy detectors" (EDI and ED2). 
The mass analyzer consists of a second cylindrical 
electrostatic deflection system suspended in a 978 G magnetic 
field. The ions that exit this region, through the image slit, 
are detected by a high-current electron multiplier (Shelley et 
ai., 1985). This instnunent was active from August 26, 1984 
until it failed on April 4, 1985. 

The EBEM consisted of eight independen 18 oo permanent 
magnet electron spectrometers. Electrons entered the EBEM 
through a 5° full angle collimator and were then deflected 
through 180° by a permanent magnet. Figure 20. They were then 
focused onto an exit aperture, that defined the allowed 
momentum range, and were then detected by a channel electron 
multiplier. (Shelley et aJ., 1985) 

Both the ion and electron data for the AMPTE/CCE HPCE were 
processed into pool files. These were the data files used for 
survey. These pool files consist of data arranged in 6.5 
minute bins from 0000 to 2400 universal time. There is a 
separate data file for each day's data and each file contains 
both electron and ion data for that day. (Shelley et al., 
1985) 
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For this work, the ion flux measurements from the energy 
detectors (ED) are used, since we were not interested in 
differentiating between and He* for this survey. The pool 
data are sorted, by time, into 18 logarithmically spaced 
energy bins. The lowest four channels use the "RPA” mode data. 
The bulk of the data which are availcdsle in the pool file have 
only the fourth RPA channel, which provides an integral 
measurement from approximately 30 to 150 eV, with a weighted 
center at 50 to 65 eV. The remaining channels extend up to 
17keV/e. The ion flux is also sorted by time verses pitch 
angle, with pitch angle bins from 0° to 90® in increments of 
lOo. (Shelley et aJ., 1985) 

The electron data is likewise sorted into 8 energy bins 
from 50 eV to 25 keV and by pitch angle from 0° to 90® in 10® 
increments, each also verses time. The energy channels for the 
ion electrostatic energy analyzer and for the EBEM are given 
in TABLE I and II. Data was also collected for ions species 
verses time verses energy and for ion species verses time 
verses pitch angle (Shelley et aJ., 1985). This aspect of the 
data was not used in this paper. 
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TABLE I. ION ENERGY CHANNELS IN THE HPCE ON THE AMPTE/CCE 


Energy Channel 

1 

2 

3 

4 


5 

6 

7 

8 

9 

10 
11 
12 

13 

14 

15 

16 

17 

18 


Energy of channel 
Center (keV/e) 

0.0014 

0.0067 

0.020 

0.050 

0.052 

0.055 

0.065 

0.240 

0.442 

0.657 

0.885 

1.127 

1.660 

2.261 

2.941 

3.709 

5.058 

6.668 

9.339 

12.75 

17.11 


Full Energy width 
(keV/e) 

0.0028 

0.0076 

0.020 

0.125 RPAl 
0.122 RPA2 
0.115 RPA3 
0.095 RPA4 
0.216 
0.209 
0.222 
0.235 
0.338 
0.567 
0.641 
0.724 
1.056 
1.480 
2.143 
3.041 
3.884 
2.600 
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TABLE II. ELECTRON ENERGY CHANNELS IN THE HPCE ON THE 
AMPTE/CCE 


Electron Detector 

CMEA 

CMEB 

CMEC 

CMED 

CMEE 

CHEF 

CMEG 

CMEH 


Energy Channel 
Center (hev/e) 

0.067 

0.150 

0.340 

0.770 

1.74 

3.94 

8.89 

20.1 


Full Energy Width 
(kev/e) 

0.051 

0.126 

0.284 

0.643 

I. 45 
3.28 
7.41 

II. 7 












III. OBSBRVlkTZOHS 


A. STATISTICIO. SURVEY 
1. Data Analysis 

Data collected from August 26, 1984 until December 6, 
1985 were processed and analyzed. However, ion data were only 
available through April 4, 1985. This was due to the failure 
of the ion-mass spectrometer on that day. The stop date for 
the analysis was chosen because the satellite completed one 
precession around the Earth, starting from August 26, 1984 and 
ending on December 6, 1985. 

For each day of data, two different types of 
spectrograms were produced for both the ions and electron 
data. The first type was an energy channel-time spectrogram, 
for pitch angles in the 80-90 degree range. The second was a 
pitch angle-time spectrogram. The energy channels used were 
the ion fifth and sixth channels (240 and 442 eV) and the 
electron third and fourth channels (340 and 770 eV). 

These spectrograms were examined for periods when the 
satellite encountered the magnetopause. The spectrograuns were 
also inspected to find periods of very "choppy” data and 
periods when the instriiment was undergoing diagnostic testing. 
Data in these periods were then removed from the data to 
ensure that the analysis would not be contaminated by such 
artifacts. 
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These edited data files were then surveyed to produce 
local time vs Mcllvain L probability distribution plots. The 
results for ions and electrons were plotted separately then 
compared. 

2. Local Time - Mcllvain L Surveys 
a. ICN SURVEY 

Figure 21 shows a plot of the probability 
distribution for eguatorially trapped ions using the lowest 
criteria from the previous chapter, namely the ion flux in the 
80°“90° pitch angle bin, at 240 eV, be greater than 10‘ 
ions/(cm^s sr), that the anisotropy be greater than 1.5, and 
that the ions be observed within 10“ latitude from the 
magnetic equator. 

The grey scale for the results runs from 0% to 60% 
with zero being white and 60% being black. The coverage plot's 
grey scale ranges from 0 to 200 counts and from white to black 
respectively. The scales are allowed to saturate (peak 
coverage was approximately 600 samples at apogee). Note that 
the 1800 to 2400 local time sector's zero occurrence 
probability is due to lack of coverage. 

The contents of Figure 21 are alternately presented 
as a surface plot in Figure 22. This plot has x and y axes of 
local time and L. The z axis is probability of occurrence, from 
0 to 100%. A contour plot is also displayed as part of this 
figure to facilitate the reading of the surface plot. Again, 
the 1800 to 2400 local time sector was not sampled. 
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Figure 21. Trapped Ions (240 eV) - Flux gt 10% Anisotropy gt 
1 .5, Maglat gt 10° 
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Figure 22. Trapped Ions (240 eV) - Flux gt 10% Anisotropy g 
1.5, Maglat gt 10° - Surface Plot 






The high probability (greater than 30%) region 
these ions starts at 0500 local time at an L of 3.5. As local 
time increases so does the L value of the peak probability. At 
1400 local time the maximum L value, 8, is reached. At this 
point the probability appears to drop off sharply in L as 
local time moves toward dusk. This, however, may be an 
artifact brought about by the low coverage after 1800 local 
time. 

In Figures 23 and 24 the results of raising the 
selection criteria for a trapped ion distribution are shown. 
In Figure 23 the selection criteria are flux greater than 10% 
anisotropy greater than 1.5, and measurements within 5° of the 
magnetic equator. In Figure 24 the selection criteria are flux 
greater than 5 x 10% anisotropy greater than 2.0, and 
measurements within 5° of the magnetic err’ator. The local time 
dependence is similar in these three figures. 

Figure 25 shows a plot of the probability 
distribution for equatorially trapped ions for the sixth 
energy channel (i.e. 442 eV). Here the selection criteria is 
a flux greater than 10* ions/(cm* s sr), an anisotropy greater 
than 1.5, and a magnetic latitude between within 5“ of the 
magnetic equator. The grey scales are the same as in the 
previous two figures. Figure 26 is a surface plot of the same 
data. 

The high probability (greater than 30%) region for 
these ions starts at 0500 local time at an L of 3.5. Again, as 
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Figure 23. Trapped Ions (240 eV) - Flux gt 10% Anisotropy gt 
1.5, Maglat gt 5° 
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Figure 24. Trapped Ions (240 eV) - Flux gt 5x10% Anisotropy 
gt 2.0, Maglat gt 5° 


57 






















AM PTE SURVEY 

IONS FROM ESA 442 EV 

Anisotropy gt 1.5 Flux gt 1.00E+03 MogLot It 5.0 


PROBABILITY DISTRIBUTION 



0 2 4 6 8 10 12 14 1618 20 22 24 
Local time 


Figure 25. Trapped Ions (442 eV) - Flux gt 10% Anisotropy gt 
1.5, Maglat gt 5° 


58 


















as local tiime increases so does the L value of the peak 
prob 2 d>ility. At 1400 local time the maximum L value, 8, is 
reached and then the probability drops off shaxply, probedsly 
due to the lack of data after 1800. 

The survey of the additional energy channels for 
the ions did not dramatically change the results found by 
Braccio. Figures 27 and 28 are the probability distribution 
plot and surface plot from Braccio (1991) for ions in the 
fourth RPA energy channel (30 - 150 eV) with the criteria of 
flux in the 80°-90° pitch angle bin be greater than 10® 
ions/(cm^s sr), that the anisotropy be greater than 1.5, and 
that the measurements be within 10° latitude of the magnetic 
equator. It is found that the high probability region (greater 
than 45%) for the ions started at 0500 local time and in the 
area of L = 3.5. His maximum L value is on average about 8.0 
and occurs generally around 1400 local time. At the maximum 
the probability appears to drop off sharply in L as local time 
moves toward dusk, similar to the results above. 

The apparent decrease in probability of occurrence 
for ions in the region of 1200 to 1400 local time found by 
Braccio is readily evident in both energy channels surveyed, 
which seems to dampen Braccio's hypothesis that the higher 
energy ions inhabited this region. 
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Figure 27. Trapped Ions (30-150 eV) - Flux gt 10^, Anisotropy 
gt 1.5, Maglat gt 10® 
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Figure 28. Trapped Ions (30-150 eV) - Flux gt 10*^, Anisotropy 
gt 1.5, Maglat gt 10® - Surface Plot 
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b. ELECTRON SURVEY 


Figure 29 shows the probability distribution for 
trapped electrons with an energy of 340 eV meeting the 
criteria of flux greater than 5 x 10* electrons/(cm* s sr), 
anisotropy greater than 1.5, and measurements within 10° of 
the magnetic equator. Figure 30 shows the seune type of plot 
for 770 eV electrons meeting the same criteria. Figure 31 
presents this 770 eV data as a surface plot. It is readily 
apparent that the electron distribution is vastly different 
from that of the ions. There is no obvious L versus local time 
trend for the electron probability distribution, similar to 
the results found by Braccio, who also fowmd the electron 
probability distribution to be conical in shape with the 
region of highest probability being between 1000 and 1100 
local time and for L values around 6.5 - 7.0, as shown in the 
probability and surface plots. Figures 32 and 33, taken from 
Braccio (1991). Here, the criteria is flux greater than 5 x 
10* electrons/(cm* s sr) , anisotropy greater than 1.5, and 
measurements within 10° of the magnetic equator. 

The cone is not completely symmetrical. The 
probability of occurrence increases to its peak value more 
gradually, as a function of local time, than it decreases. 
This characteristic does not seem to alter even when changing 
the selection criteria or energy channels and is the same as 
that found by Braccio. The statistical analysis of the 
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Figure 29. Trapped Electrons (340 eV) - Flux gt 5x10* 


Anisotropy gt 1.5, Maglat gt 10 
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Figure 30. Trapped Electrons (770 eV) 
Anisotropy gt 1.5, Maglat gt 10° 
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Figure 32. Trapped Electrons (150 eV) - Flux gt 5x10' 
Anisotropy gt 1.5, Maglat gt 10® 
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Figure 33. Trapped Electrons (150 eV) - Flux gt 5x10 
Anisotropy gt 1.5, Maglat gt 10® - Surface Plot 
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coarsely gridded survey files is now extended using more 
detailed looks at the raw data files. 

B. DETAILED ANALYSIS CASE STUDIES 
1. Data Analysis 

Detailed analysis of the electron data was conducted, 
using data at 150 eV and 340 eV. ion data were not readily 
available. Data were collected for the days in TABLE III. Only 
days 84/336 and 85/002 will be presented in the thesis. The 
days listed in TABLE III were selected from the set of days 
where ••large” electron events occurred, as found in Braccio 

Table III. DETAILED ANALYSIS DATA COLLECTION DAYS 


Day 

Local Time 
Range 


D.t 

84/266 

0941-1634 

26 

-21 

84/283 

0510-1538 

26 

-24 

84/296 

0755-1420 

39 

-55 

84/315 

0655-1349 

23 

-11 

84/336 

0533-1249 

22 

-24 

85/002 

0359-1418 

24 

-3.2 

85/039 

0129-0834 

28 

-24 

85/227 

1440-1938 

20 

-21 

85/238 

1330-2330 

22 

-10 

85/264 

2002-2156 

29 

-19 

85/301 

1618-2049 

7 

-17 

85/320 

1116-1837 

19 

-18 

85/337 

1021-1742 

17 

-31 
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(1991) and extended here, to obtain reasonably complete local 
time coverage. Special emphasis was given to the morning 
sector (0600 to 1200 local time) . Each segment of data 
considered covered the period between plasmapause encoiinters, 
roughly 15 hours of data per day. Also, listed with the days 
and times in TABLE III are the magnetic activity indices, EKp 
and Djt, for the respective time period. 

Data were initially considered in one hour increments. 
For each hour and each of the energy channels surveyed, 150 
and 340 eV, a pitch angle vs UT time spectrogram was 
generated. Figure 34 is a compilation of the one hour pitch 
angle vs time spectrograms in reduced form for 150 eV 
electrons on day 84/336. The 150 eV spectrograms for the 
additional days surveyed can be found in the Appendix. The 
grey scale for the spectrogram runs from 7.2 to 9.0 on a log 
scale of the flux. Pitch angle runs from 0' - 180° and UT from 
0000 to 1400. Also, plotted on the x-axis is the corresponding 
local time (LT), Mcllwain L value and magnetic latitude. Note 
the high flux at 90 degrees from 1000 to 1400 UT. Note how the 
equatorially trapped distribution fades as the satellite 
approaches perigee (the plasmapause at L " 3.9.) This 
characteristic can be seen in a number of the survey 
spectrograms as shown in Appendix. The data were analyzed by 
considering the log of distribution function as a function of 
cos^o. Figure 35. This functional form should produce straight 
line segments, if the data are well represented as a 
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Figure 34. Pitch Angle verses Time Spectrogram - Day 84/336 
Energy 150 eV 
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Figure 35. Electron Pitch Angle Distribution 
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bi-Maxwellian (Scott, 1991.) From these segments values for 
Tperp and Tpari Can be estimated, and various fliix ratios 
obtained. Much of this work is based on consideration of how 
the data (log f) vary with pitch angle (cos^a) . 

From this figure a distinct break in the vicinity of 
60 degrees can be seen. At lower pitch angles the isotropic 
background dominates the pitch angle distribution. Because of 
this effect the data analysis range was limited to between 60° 
and 90°, a restriction that more accurately considers the 
phenomenon of interest. This restriction, however, has the 
effect of limiting the data set which can be analyzed. These 
equatorial pitch angles (60° to 90°) do not map far from the 
equator. A restriction of +-6° magnetic latitude was chosen, 
which allows a suitable range of (equatorial) pitch angles, 
and allows a reasonable range of events. 

Fluxes were obtained in three ways; 

1. Average Flux (60° - 70°): the flux in the 60° to 70° pitch 
angle bin at the location of the satellite. 

2. Average Flux (80° - 90°) ; the flux in the 80° to 90° pitch 
angle bin at the location of the satellite. 

3. Fitted Flux (60° - 90°); the fitted flux was derived by 
taking a least squares fit of the data at pitch angles 
between 60° and 90° on a plot of the log of the distribution 
function versus cos^a. This fit was used to project the 
values which would have been obtained at the magnetic 
equator, effectively assuming a bi-Maxwellian form. 
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Two different anisotropies were calculated and are 
defined as follows: 

1. Average Anisotropy: the ratio of the average fluxes in 
the 80° to 90° pitch angle bin to those in the 60° to 70° 
pitch angle bin. 

2. Equator Anisotropy: the fluxes at the latitude of 
measurement mapped to the equator and a ratio taken for the 
flvix at 90° versus the flux at 69°. 


Two temperatures were calculated for each energy 

channel: 

i* *^perp • the characteristic temperature in the 
perpendicular direction, with respect to the magnetic field 
line, of a bi-Maxwellian distribution. Calculated from the 
difference in the log of the distribution functions for 150 
eV and 340 eV at a pitch angle of 0°. The calculated value 
is the same for both energy channels. 

2. Tpgji : the characteristic temperature in the parallel 
direction, with respect to the magnetic field line, of a 
bi-Maxwellian distribution. This is based on the calculated 
Tperp and the slope of the pitch angle distribution at the 
respective energy. The calculated value is different for 
each energy channel, reflecting the non-Maxwellian element 
of the energy distribution. 


Both temperatures and the density were calculated from 
a bi-Maxwellian form for the distribution function. This 
produces valid and well behaved results if the distribution 
can be characterized as a bi-Maxwellian. If it is not well 
described as a bi-Maxwellian the results will be invalid and 
not well behaved. 

The ratio of Tp^rp to Tpa^-j was also calculated for the 
data. This ratio was used to further characterize a trapped 
distribution. For the purposes of this paper the ratio had to 






exceed three for a distribution to be characterized as 
eguatorially trapped. This corresponds to a density decrease 
factor of 20% at 10° magnetic latitude in the aibsence of 
parallel electric fields. (Scott, 1991) 

2. Day 84/336 Case Study Trapped Distributions 
a. 150 eV 

Figure 36 shows the plots of the previously defined 
fluxes vs UT time at 150 eV. The solid horizontal line plotted 
with the fluxes in the bottom two plots is for the flux value 
of 5x10* (cm^ s sr)*^ , the minimum threshold established in the 
last chapter. Also, plotted on the x-axis are the local time, 
Hcllwain L and the magnetic latitude. 

Figure 37 shows the plots of the previously defined 
anisotropies versus time at 150 eV. The solid horizontal line 
plotted on each is the anisotropy threshold of 2.0, defined in 
the previous chapter as the minimum anisotropy allowed to 
characterize a trapped distribution. Note the periods that the 
anisotropies exceed the threshold correspond to the periods of 
increased flux. Also, shown in Figure 37 is the plots of Tperp 
vs time, Tpa^i vs time, and the ratio Tperp to Tpan vs time. The 
solid horizontal line plotted on the plots for temperature 
ratios vs time is the ratio threshold level of three defined 
in the last chapter. At the time the temperature ratio exceeds 
three perpendicular temperature is on average 70 eV while the 
parallel temperature drops from values of 20 eV to 10 eV. The 
peak of the ratio has a value of 6.0. For this day a 
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Figure 36. Fluxes verses Time - Day 84/336, Energy 150 eV 
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eguatorially trapped electron distribution appears to occur 
for a Mcllwain L value of about 6 to 7.5. Note that the lower 
anisotropy, 1.5 (the dashed line in Figure 37), used in the 
statistical survey would give a wider L extent, from L = 3 to 
9. With this lower anisotropy and the corresponding L range as 
a boundary, for where eguatorially trapped electron 
distributions exist, a flux value of 5x10^ and a temperature 
ratio of two can be interpreted as new criterion thresholds 
for eguatorially trapped electrons, as seen by the dashed 
lines in the flux and temperature ratio plots of Figures 36 
and 37. 

Figure 38 shows the plot of the average anisotropy 
versus the eguator anisotropy and the plot of temperature 
ratio versus eguator anisotropy. The average anisotropy is 
slightly higher than the fitted eguatorial anisotropy as can 
be seen 'when compared to a line of slope one (the line on the 
plot.) In the plot of temperature ratio vs eguator anisotropy 
the data starts to exceed the ratio criterion threshold of 
three at an eguator anisotropy of about 2.0, if the high end 
of the data points is taken. This tends to support our 
definition for an eguatorially trapped electron distribution 
using flux ratios, and shows that a bi-Maxwellian description 
gives comparable results. 

Using the criteria for eguatorially trapped 
electron distributions defined in the previous chapter of a 
minimum flux for 150 eV electrons of 5x10® (cm^ s sr)S an 
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Figure 38. Average Anisotropy and Temperature Ratio verses 
Equator Anisotropy - Day 84/336, Energy 150 eV 







average anisotropy of 2.0, a measurement within 10° of the 
magnetic equator, and the additional criterion for a ratio of 
Tperp vs Tp„i to be greater than three an equatorially trapped 
plasma can be seen to start at around 0942 local time at a L 
value of 8.7 and at magnetic latitude of -2.4°. As L decrease 
with time all anisotropies increase to a peak at around 1050 
local time and a L value of 6.9 and a temperature ratio of 
6.0. All anisotropies decrease with time after that until the 
trapped distribution ends at around 1147, L of 5.30 and 
magnetic latitude of 2.4°. 
b. 340 eV 

Figure 39 shows the plots of the previously defined 
fluxes vs UT time at 340 eV. The solid horizontal line plotted 
with the fluxes in the bottom two plots is for the flux value 
of 5x10* (cm^ s sr)'^ , the minimxim threshold established in the 
previous chapter. Also, plotted on the x-axis are the local 
time, Mcllwain L and the magnetic latitude. 

Figure 40 shows the plots of the previously defined 
anisotropies versus time at 340 eV. The solid horizontal line 
plotted on each is the anisotropy threshold of 4.4, defined in 
the previous chapter as the minimum anisotropy allowed to 
characterize a trapped distribution. Also, shown in Figure 40 
is the plots of Tpe,,p vs time, Tpari vs time, and the ratio Tpe^p 
to Tpari vs time. The solid horizontal line plotted on the 
plots for temperature ratios vs time is the ratio threshold 
level of three defined in the last chapter. At the time the 
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Figure 39. Fluxes verses Time - Day 84/336, Energy 340 eV 
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Figure 40. Anisotropies, Temperature Ratio and Temperatures 
verses Time - Day 84/336, Energy 340 eV 








temperature ratio exceeds three perpendicular temperature is 
on average 70 eV while the parallel temperature drops to 
values of 23 eV. At this energy a trapped electron 
distribution appears to occur for a more limited Mcllwain L 
range, at the higher end of the L range defined by the 150 eV 
data. 

Using the criteria for trapped electron 
distributions defined in the previous chapter of a minimum 
flux for 340 eV electrons of 5x10® (cm^ s sr)'S an average 
anisotropy of 4.4, a measurement within 10° of the magnetic 
equator, and the additional criterion for a ratio of Tp^^p vs 
Tpari to be greater than three an equatorially trapped plasma 
can be seen to start at around 1031 local time at a L value of 
7.5 and at a magnetic latitude of 0.8°. The trapped 
distribution does not extend very far ending at 1041 at L 
equal to 7.2 and magnetic latitude of 1.2°. A more realistic 
anisotropy ratio of " 2 and a temperature ratio of 2 , as 
shown by the dashed line in both plots of Figure 40 would give 
a latitude extent more in line with respect to the 150 eV 
results. 

3. Day 85/002 Case study Trapped Distributions 
a. 150 eV 

Figure 41 shows the plots of the 150 eV fluxes vs 
time for 85/002. The solid horizontal line plotted with the 
fluxes in the bottom two plots is for the flux value of 5x10® 
(cm^ s sr)''. Also, plotted on the x-axis are the local time. 
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Mcllwain L and the magnetic latitude. Comparison with the 
anisotropy suggests a flux level of 5x10^ would be a more 
accurate indicator on this day, as shown by the dashed line in 
the plot. 

Figure 42 shows the plots of the previously defined 
anisotropies versus time at 150 eV. The solid horizontal line 
plotted on each is the anisotropy threshold of 2.0, defined in 
the previous chapter as the minimum anisotropy allowed to 
characterize an eguatorially trapped distribution. Note the 
periods that the anisotropies exceed the threshold 
corresponding to the periods of increased flux (> 5x10’) . 
Also, shown in Figure 42 are the plots of Tperp vs time, Tpari vs 
time, and the ratio Tperp to Tpari vs time. The solid horizontal 
line plotted on the plots for temperature ratios vs time is 
the ratio threshold level of three. At the time the 
temperature ratio exceeds three perpendicular temperature is 
on average 80 eV while the parallel temperature drops to 
values of 20 eV. For this day a trapped electron distribution 
appears to occur for a Mcllwain L value of about 5 to 7. Note 
the lower anisotropy ratio of 1.5 (the dashed line in Figure 
42), used in the statistical survey would give a wider L 
extent, from L = 5 to 8. With this lower anisotropy and the 
corresponding L range as a boundary, for where eguatorially 
trapped electron distributions exist, a flux value of 5x10’ 
and a temperature ratio of two can be interpreted as new 
criterion thresholds for eguatorially trapped electrons. 


84 




ur 

6.0 

n 

8.0 

1 

10:0 

1 

1Z0 

r 

14i) 

LT 

7J 

8.0 

a7 

u 

14.1 

L 

U 

&7 

7.7 

5.7 

ZO 

X -4.4 

-1i) 

zo 

i6 

U 


Figure 41. Fluxes verses Time - Day 85/002, Energy 150 eV 
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as seen by the dashed lines in the flux and temperature ratio 
plots of Figures 41 and 42. 

Using the criteria for trapped electron 
distributions defined in the previous chapter of a minimum 
flux for 150 eV electrons of 5xl0‘ (cm^ s sr)'‘, an average 
anisotropy of 2.0, a measurement within 10° of the magnetic 
equator, and the additional criterion for a ratio of Tp^rp vs 
Tpari to be greater than three a equatorially trapped plasma 
can be seen to start at around 0900 local time at a L value of 
7.2 and at magnetic latitude of 2.6°. As L decreases with 
time the anisotropies increase to a peak at around 0909 local 
time and an L value of 6.9. The anisotropies decrease with 
time after that until the trapped distribution ends at around 
1003, L of 5.30 and magnetic latitude of 3.6“. 
b. 340 eV 

Figure 43 shows the plots of the fluxes vs time at 
340 eV. The solid horizontal line plotted with the flxixes in 
the bottom two plots is for the flux value of 5x10® (cm* s 
sr)*. Also, plotted on the x-axis are the local time, Mcllwain 
L and the magnetic latitude. 

Figure 44 shows the plots of the anisotropies 
versus time at 340 eV. The solid horizontal line plotted on 
each is the anisotropy threshold of 4.4. Also, shown in 
Figure 44 are the plots of Tp„p vs time, Tp^^i vs time, and 
theratio Tp^rp to Tpa^i vs time. For the short interval where the 
temperature ratio exceeds three the perpendicular temperature 
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Figure 43. Fluxes verses Time - Day 85/002, Energy 340 eV 


88 






UT 

— 1- 

-1- 

to 

-1- 

10.0 

-1- 

12.0 

- 

14^ 

LT 

7J 

to 

&7 

S.8 

14.1 

L 


B.7 

7.7 

5.7 

2.0 

X 

-4.4 

-^0 

20 

3.6 

1.2 


Figure 44. Anisotropies, Temperature Ratio and Temperatures 
verses Time - Day 85/002, Energy 340 eV 
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is on average 90 eV while the parallel temperature drops to 
values of 30 eV. Again, the eguatorially trapped distribution 
is much more limited in its appearance at 340 eV. Indeed the 
criterion of a ratio of 4.4 is not met. However, if we lower 
the anisotropy criterion to 3.0 (the dashed line in the 
anisotropy plots of Figure 44) a trapped distribution would be 
defined to occur around a Mcllwain L value range of 6 to 7.5 
at around 0901 local time at a L value of 7.1 and at magnetic 
latitude of 2.8°. The trapped distribution does not extend 
very far ending at 0928 at L equal to 6.3 and magnetic 
latitude of 3.4°. This corresponds to a temperature ratio ' 2 
(the dashed line in the temperature ratio plot of Figure 44) 
4 . Detailed Analysis Trapped versus Non-Trapped 

Distributions 

In Figure 34 there was evidence of two areas of 
trapped electron distributions for day 84/336. The first was 
presented using Figures 36 through 38, which was an 
eguatorially trapped distribution as defined by our criteria. 
The other area that presents itself in Figure 34 is the time 
period between 0553 and 0700 local time. The broader 
distribution seen earlier in the day fail our criteria. Figure 
45 and 46 show plots like previously of the fluxes, 
anisotropy, temperatures and temperature ratios, all verses 
time for the entire orbit on 84/336 at 150 eV. Though both 
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Figure 45. Fluxes verses Time - Day 84/336, Energy 150 eV 
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Figure 46. Anisotropies, Temperature Ratio and Temperatures 
verses Time - Day 84/336, Energy 150 eV (Entire Orbit) 
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time periods show similar increases in flux the anisotropies 
for the second time period barely exceeds 1.3 and the 
temperature ratio never rises above 2.0. 

Figure 47 is a plot of the logiof versus pitch angle 
for 150 eV electrons from different days and hours in the 
detailed analysis. The plot is organized from bottom to top to 
show a progression of non-trapped to trapped electron 
distributions. Plot a and b show the pitch angle distribution 
for periods when no trapped electron distribution is present. 
Plot c and d show the beginnings of trapped electron 
distributions. Plots e and f show the narrowest equatorially 
trapped electron distributions observed in the detailed 
analysis. Note that as the trapped distribution gets stronger 
the peak at 90° gets progressively more pronounced, indicating 
that the electrons are equatorially trapped. This 
characteristic suggests that the data may be fitted with a bi- 
Maxwellian distribution function as was done. This fit 
performed well and it was quite evident it did not fit well 
indicating that the distribution was not trapped. This 
criteria allows us to distinguish between (simply) trapped 
distributions (e.g., loss c<nes), versus equatorially trapped 
distributions. 
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Figure 47. Flux verses Pitch Angle - Various Days 













ZV. DISCnSSZOM 


A. STATISTICAL SURVEY 

The probability distributions for equatorially trapped 
ions and electrons show clear localization for regions of high 
occurrence probability (>= 50%). This is particularly true of 
electrons. The electron distribution has a very localized peak 
at 1000 “ 1100 local time, L*= 6.0 - 6.5, and within 5 degrees 
of the magnetic equator. This is seen regardless of the 
selection criteria used in the surveys. This agrees with the 
results found by Braccio. The survey also showed the 
ec[uatorially trapped electrons were primarily located outside 
the traditional plasmasphere, and were primarily confined to 
the dawn to noon sector. The peak in the electron distribution 
occurred near geosynchronous orbit. 

The results for trapped ion distributions, also, agree 
closely with Braccio's results with the exception of his 
explanation for the decrease in probability between 1200 - 
1400. He suggested that the decrease in probability is due to 
the average energy of the equatorially trapped ions going 
above the 150 eV energy he surveyed, in this local time 
region. Both the 240 eV and 442 eV ions surveyed showed the 
same decrease in ion probability distribution for this local 
time region. It maybe possible that Braccio's hypothesis is 
correct. Sagawa et al (1987) did not show a decrease in 
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probability for the region from 1200 to 1500 local time. 
However, Sagava et al. used summary results from approximately 
10 eV to 1 keV, whereas Braccio's survey and this survey 
covered 50 to 442 eV. To accurately determine if higher energy 
ions inhabit this region of decreased probability the 
additional energy channels of 657, 885 and 1127 eV ought to be 
surveyed. The ion survey reinforces Braccios* and other 
previous work with ISEE-1 and DE-1, with peak ion flux inside 
the plasmapause. 

B. DETAILED ANALYSIS CASE STUDIES 

In the analysis only 13 days were looked at, therefore a 
statistical analysis isn't appropriate. However, the orbital 
coverage. Figure 48, of the days analyzed was widespread 
enough to see patterns or tendencies emerge in the data. 

Figure 49 is a plot of where the 150 eV data from the 13 
days that met the criteria to be considered a trapped electron 
distribution occurred; namely a flux greater than 5.0x10* (cm^ 
s sr)"^, an anisotropy greater than 2.0, a temperature ratio 
greater than 3.0 and within 10® of the magnetic equator. The 
data is plotted in the x-y plane of the AMPTE satellite orbit. 
The trapped distributions start at dawn at a mean L value of 
about 3.8 and increase in orbital density and L to a maximum 
between 0900 and 1100 local time and an average L of 7.0. The 
distributions then drop off rather rapidly with only smaller 
occurences of equatorially trapped distributions in 
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mid-aftemoon at L values in excess of 7.0 and at dusk when 
the distributions move back into the smaller L values of 5.0 
to 6.5. 

Braccio (1991) in his statistical survey found a similar 
general pattern to exist for trapped electron distributions at 
the 150 eV energy range. He found that the electron 
distribution had a very localized peak at 0900 local time, L 
= 6, and within 5° of the magnetic equator. The difference in 
the results can be attributed to his survey considering all 
days of operation the AMPTE satellite while this thesis only 
considered a small portion of the operational life. 

The tendency for the trapped electron distributions to 
move out further in L value starting at dawn until late 
afternoon and then move back in near dusk can be explained in 
terms of plasmapause location. Figure 50 is taken from 
Williams et al. (1981). It shows the location of the 
plasmapause as a function of local time. The shape of this 
plot mimics the behavior of the trapped electron distributions 
in the case studies. Then as the plasmapause moves out during 
the day so does the location of the trapped electrons. 

Braccio (1991) found that trapped electrons seem to be 
excluded from regions of high trapped ion probabilities and 
vice versa. He, also, found that the trapped electron 
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distributions probabilities drop off by local noon at the 
beginning of the peak of the trapped ion probabilities, which 
he found to peak between 1400 and 1500 local tine at L = 8. 
Horowitz et ai. (1981) observed that trapped ion distributions 
tended to occur along the plasnapause boundary in the dawn to 
dusk region with the majority just inside the boundary. 
Therefore Braccio postulated that if trapped ions occur at the 
plasnapause, then as the altitude of the plasnapause rises, in 
the course of a day, the region where trapped electron 
distributions usually occur is saturated with trapped ions. 
This ion saturation seems to, in turn, lead to a disruption in 
the process by which the trapped electrons are produced. This 
approach could be used to explain the gap in trapped electron 
distributions found in the case studies. Figure 49, between 
1200 and 1700 local time inside of L » 7. 

Figure 51 is a plot of the occurrence of 340 eV data in 
the 13 days that met the criteria to be considered a trapped 
electron distribution; namely a flux greater than 5.0x10® (cm^ 
s sr)‘^, an anisotropy greater than 4.4, a temperature ratio 
greater than 3.0 and within 10° of the magnetic equator. The 
data is plotted in the x-y plane of the AMPTE satellite orbit. 
There appears to be no major grouping of trapped electron 
distributions at this energy level. This is probably the 
result of the higher anisotropy requirement of 4.4. The 
trapped distributions are localized to 0900 to 1100 local time 
at L of 7.0 and at dusk from 1700 to 1800 local time at L=5.0. 
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A comparison of the two different anisotropies used to 
characterize the data found that they tracked reasonedsly well 
with the equator anisotropy being on average 1.075 times the 
equator anisotropy in the dawn to noon region and only 0.59 
times the equator anisotropy in the noon to post dusk region. 
The difference in anisotropy ratios for the equator anisotropy 
between the morning and the afternoon local time could be due 
to a lack of samples taken in the afternoon region that were 
characterized as trapped. 

Over 28,000 samples taken in 13 orbits are summarized in 
Figure 52, the plot of Tp^,ip/Tpaii vs Equator Anisotropy for 150 
eV electrons. This figure shows a different technique for 
defining anisotropy of a equatorially trapped electron 
distribution and suggests that there is a monotonic 
relationship between a temperature ratio greater than 3 and an 
anisotropy of 2.0 for 150 eV and 4.42 for 340 eV. For all days 
the threshold of 3 for a temperature ratio begins to be 
exceeded at the anisotropies of about 2.0 for 150 eV and about 
4.42 for 340 eV. 
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Figure 52. Temperature Ratio verses Equator Anisotropy - All 
13 Detailed Analysis Case Study Days, Energy 150 eV 
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V. COHCLnSIOMS 


A. STATISTICAL ANALYSIS 

The location of the eguatorially trapped plasmas are 
species dependent. The ions and electrons show a different 
local time dependence in the location of their occurrence 
probability peak. Electrons show a uniform high probability 
distribution centered at 0900 local time and an L value of 6. 
The fact that trapped electrons begin to be seen at local dawn 
lead to speculation that their existence is dependent on 
photoelectron emission form the Earth's ionosphere. The shape 
of this distribution is basically conical, however, it drops 
off more rapidly than it increases, with respect to local 
time. 

Ions, meanwhile, have a probability distribution that 
shows a strong L dependence on local time. The high trapped 
probability region begins at local dawn, for L approximately 
4, and rises to a maximxim at between 1400 and 1500 local time 
with an L value of 8. The distribution then drops quickly in 
altitude as local time increases. Additionally, there is a 
region of decreased probability in the afternoon sector for 
this data that was speculated by Braccio (1991) to be 
inhabited by higher energy ions, but still is present for the 
higher energies surveyed. 

The overall shape of the high trapped ion probability 
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region mirrors that of the location of the plasmapause. This 
suggests that the trapped ion distribution is linked to the 
plasmapause. This would confirm Horwitz et al. 's (1981) 
observations that 'pancake' distributions often occur in the 
vicinity of the plasmapause. 

Additionally, it has been observed that trapped electrons 
seem to be excluded from regions of high trapped ion 
probabilities and vice versa. If the trapped ions actually do 
occur at the plasmapause, then as the altitude of the 
plasmapause rises, in the course of a day, this exclusion 
would effectively create a barrier that restricts the trapped 
electron distribution to the dawn to noon sector. 

B. DETAILED ANALYSIS CASE STUDIES 

The results of the detailed analysis case studies in 
general agree with the above results found by Braccio (1991). 
Trapped electron distributions are concentrated in the dawn to 
noon ragion with a peak around 1000 at an L value of 7. What 
trapped distributions occur in the afternoon are generally at 
higher L values until late afternoon when they move back into 
L values around 5 to 6 at 1800. This follows with the results 
found by Braccio (1991) who found that large trapped electron 
distributions are exclude from areas of high trapped ion 
distributions that peak in the 1400 to 1500 local time region. 
This also suggests a connection between the plasmapause and 
trapped electron distributions. 
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The threshold flux level used in the statistical survey 
was too low for 150 eV electrons. A higher flux level of 5x10' 
(cm^ s sr)'^ could be used for this energy level. The use of an 
anisotropy of 1.5 and a temperature ratio of 2.0 produced 
results more in keeping with the results of the statistical 
survey. The threshold flux level for 340 eV appeared to be 
correct. 

The use of an average anisotropy at a latitude within 6° 
of the magnetic equator vice a calculated equator anisotropy 
in defining trapped electron distributions is generally pretty 
accurate. The average difference is only 7.5% between the two. 

It has been shown that equatorially trapped plasmas can be 
described by a bi-Maxwellian distribution function. With such 
a description the values of the perpendicular and parallel 
temperatures can be calculated and their ratio used along with 
an anisotropy ratio between the flux in the 80® to 90° pitch 
angle range with the flux in the 60° to 70° pitch angle range 
to define a electron distribution as equatorially trapped or 
not. It was found that a temperature ratio of three along with 
anisotropies of 2.0 for 150 eV electrons and 4.4 for 340 eV 
electrons could be used to define a trapped electron 
distribution along with the required flux and magnetic 
latitude requirements. 

However, the 340 eV anisotropy of 4.4 produce less 
reliable results than the 150 eV anisotropy with the 
theoretical anisotropy and temperature ratios for 340 eV being 
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too high. In this study Tp^^p is calculated from the bi- 
Maxwellian distribution function and this results in the same 
Tp,rp for both 150 eV and 340 eV. This was the only method 
available. However, the method produces a Tp^^p that is an 
underestimate of the perpendicular temperature for 340 eV 
particles. This suggests that the trapped electron 
distributions at this energy are not true bi-Maxwellians, but 
are more like a kappa function. 






APPENDIX (DETAILED ANALYSIS SPECTROGRAMS) 
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Figure 53. Pitch Angle verses Time Spectrogram - Day 84/266, 
Energy 150 eV 
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Figure 55. Pitch Angle verses Time Spectrogram - Day 84/296, 
Energy 150 eV 
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Figure 57. Pitch Angle verses Time Spectrogram - Day 85/002, 
Energy 150 eV 
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Figure 58. Pitch Angle verses Time Spectrogram - Day 85/039 
Energy 150 eV 
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Figure 59. Pitch Angle verses Time Spectrogram - Day 85/227, 
Energy 150 eV 
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Figure 60. Pitch Angle verses Time Spectrogram - Day 85/238 
Energy 150 eV 
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Figure 61. Pitch Angle verses Time Spectrogram 
Energy 150 eV 
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Figure 63. Pitch Angle verses Time Spectrogram 
Energy 150 eV 
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Figure 68. Pitch Angle verses Time Spectrogram - Day 84/315, 
Energy 340 eV 




















































Figure 70. Pitch Angle verses Time Spectrogram - Day 85/002, 
Energy 340 eV 
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AMPTE 85227 34-0 eV Electrons 










































































































































Figure 75. Pitch Angle verses Time Spectrogrcim - Day 85/301 
Energy 340 eV 


































Figure 76. Pitch Angle verses Time Spectrogram - Day 85/320 
Energy 340 eV 
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